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Abstract 

We present an efficient decoding algorithm for constant rate quantum hypergraph-product 
LDPC codes which provably corrects adversarial errors of weight Vt{^/n) for codes of length n. 
The algorithm runs in time linear in the number of qubits, which makes its performance the 
strongest to date for linear-time decoding of quantum codes. The algorithm relies on expanding 
properties, not of the quantum code’s factor graph directly, but of the factor graph of the original 
classical code it is constructed from. 
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1 Introduction 


A quantum CSS code is a particular instance of a quantum stabilizer code, and can be defined 
by two classical binary linear codes Cx and Cz in the ambient space F 2 , with the property that 
C Cz and C Cx- In other words, the classical codes Cx and Cz come together with respective 
parity-check matrices Hx and Hz such that the linear space Rx = Cj^ generated by the rows of 
IIx is orthogonal to the row space Rz = C^ of Hz, where orthogonality is with respect to the 
standard inner product. An error pattern is defined as a couple {ex,ez), where ex and ez are 
both binary vectors. The decoder is given the pair of syndromes ax = HxOx = Hze^ and 

decoding succeeds if it outputs, not necessarily the initial error pattern {ex,ez), but a couple of 
the form {ex + fx,&z + fz) where fx £ Rz and fz G Rx- See [H] for the equivalence with the 
stabilizer formalism and a detailed introduction to quantum coding. 

If efficient quantum computing is to be achieved, it will come with a strong error-correcting com¬ 
ponent, that will involve very fast decoders, probably in not much more than linear time in the 
blocklength n. The likeliest candidates for this task are quantum LDPC codes: in the CSS case, an 
LDPC code is simply a code whose above parity-check matrices Hx and Hz have row and column 
weights bounded from above by a constant. Among recent developments, the recent paper |12j has 
shown how fault tolerant quantum computation with constant multiple overhead can be obtained, 
and quantum LDPC codes are an essential component of the scheme, making them possibly even 
more appealing. 

It is natural to hope that the success of classical LDPC codes, both in terms of performance and of 
decoding efficiency, can eventually be matched in the quantum setting. This agenda involves two 
major difficulties, however. The first one is that coming up with intrinsically good constructions of 
quantum LDPC codes is in itself a challenge. In particular the random constructions that can be so 
effective in the classical case do not work at all in the quantum case. Indeed if one chooses randomly 
a sparse parity-check matrix Hx, then, precisely since this gives a good classical code, there are no 
low-weight codewords in the dual of the row-space of Hx and therefore an appropriate matrix Hz 
does not exist. A testament to this difficulty is given in the introduction of [18] by way of a list of 
proposed constructions of quantum LDPC codes from within the classical coding community that 
all yield constant minimum distances. Presently, the known constructions of families of quantum 
LDPC codes that come with constant rates and minimum distances that grow with the qubit 
length can be reduced to essentially three constructions. The first consists of quantum codes based 
on tilings of two-dimensional hyperbolic manifolds (surfaces) that generalize Kitaev’s toric code 
and originate in |8]. The minimum distance of these codes grows as logn, where n is the qubit 
length. A recent generalisation of this approach to 4-dimensional hyperbolic geometry m yields 
minimum distances that behave as where e is larger than some unknown constant and not more 
than 0.3. Finally, the construction [15] yields codes of constant rate with minimum distances that 
grow as These codes are perhaps the closest to classical LDPC codes in spirit, since they are 
constructed by taking a properly defined product of a classical LDPC code with itself. We note 
that presently all known constructions of quantum codes, even if they are allowed to have vanishing 
rate, fail to significantly break the barrier for the minimum distance and it is an intriguing 
open question as to whether there exist asymptotically good quantum LDPC codes (i.e. with 
constant rate and minimum distance linear in n). We also make the side remark that Gottesman 
[T^ requires, for the purpose of fault-tolerant quantum computation, constant rate LDPC codes 
with good minimum distance properties that should behave well under some sort of adversarial 
error setting. 
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The second difficulty in attempting to match the achievements of classical LDPC coding, is to devise 
efficient decoding algorithms. The vast majority of decoding algorithms developed for classical 
LDPC codes rely on iterative techniques whose ultimate goal is to take decisions on individual 
bits. In the quantum setting, taking a decision on an individual qubit is mostly meaningless: this 
is because a decoder who would try to recover the error vector exactly is doomed to fail, since it 
would be fooled by any error that spans only half a stabilizer vector (a row vector fx of Hx or fz 
of Hz) [l6|. The eventual error vector that one must look to output is therefore only defined up 
to addition of stabilizer vectors, so that there is no “correct value” for a single qubit of ex or ez 
which can always be just as well 0 or 1. 

Decoding quantum LDPC codes requires therefore additional elements to the classical tool kit. 
Surface codes mentioned above come with efficient decoding algorithms, however they are very 
particular to the surface code structure, which implies that the associated classical codes Cx and 
Cz are cycle codes of graphs: full decoding, which is NP-hard in general for linear codes, can 
be achieved for cycle codes of graphs in polynomial time (with the help of Edmonds’ weighted 
matching algorithm 0 ), and this strategy (which does not really qualify as a local technique) 
yields a decoding scheme for the quantum code that achieves vanishing error-probability for random 
errors. Unfortunately, this technique does not extend to other classes of LDPC codes, and in an 
adversarial setting is limited to correcting at most logn errors, since the minimum distance of 
surface codes of constant rate can never surpass a logarithm of the qubit length [5j . Very recently, 
an alternative decoding algorithm was proposed m for the 4-dimensional hyperbolic codes of 
Guth and Lubotzky that is devised to work in a probabilistic setting and for which its adversarial 
performance is unclear. The third class of constant rate quantum codes with growing minimum 
distance, namely the codes [l8], had no known decoding algorithm to go with it until the present 
paper whose prime objective is to tackle this very problem. 

In this paper we devise a decoding algorithm for the product codes |18j that runs in linear time 
and decodes an arbitrary pattern of errors of any weight up to a constant fraction of the minimum 
distance, i.e. for some constant c > 0. Our inspiration is the decoding algorithm of Sipser and 

Spielman m which applies to classical LDPC codes whose Tanner graph is an expander graph: 
recall that “the” Tanner graph (there may actually be several for the same code) is a bipartite 
graph defined on the set of rows and the set of columns of a parity-check matrix for the code and 
which puts an edge between row i and column j if the matrix contains a “1” in position (i, j). The 
quantum codes under consideration here are products (in a sense to be given precisely below) of a 
classical LDPC code C with itself, and we take the original code C to be an expander code. The 
resulting Tanner graphs of the two classical codes Cx and Cz that make up the quantum code are 
not strictly speaking expander graphs, but they retain enough expanding structure from the original 
code for a decoding algorithm to work. Arguably, this is the first time that an import from classical 
LDPC coding theory succeeds in decoding a quantum LDPC code from a non-constant number of 
errors in an adversarial setting. There are some twists to the original Sipser-Spielman decoding 
scheme however, since it guesses values of individual bits and we have pointed out that this strategy 
cannot carry through to the quantum setting. The solution is to work with generators rather than 
qubits: the generators are the row vectors of Hx and Hz (thus called because they correspond to 
generators of the stabilizer group of the code). At each iteration, the decoding algorithm looks for 
a pattern of qubits inside the support of a single generator that will decrease the syndrome weight. 

Our results also have some significance in the area of local testability. Locally Testable Codes (LTC) 
play a fundamental role in complexity theory: they have the property that code membership can 
be verified by querying only a few bits of a word m- More precisely, the number of constraints not 
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satisfied by a word should be proportional to the distance of the word from the code. Given their 
importance, it is natural to ask whether a quantum version of LTC exists, and to investigate their 
consequences for the burgeoning field of quantum Hamiltonian complexity, which studies quantum 
satisfaction problems [9]. 

Quantum LT codes were recently defined in [2], and these hypothetical objects are mainly charac¬ 
terized by their soundness (or robustness) R{5), i.e. the probability that a word at relative distance 
6 from the code violates a randomly chosen constraint. 

While we do not exhibit such quantum LT codes here, we construct codes which are robust for 
errors of reasonably low weight, up to a constant fraction of the minimum distance: see Corollary 
below. Reaching beyond the regime of low weight errors appears to be much harder since it 
is well-known that the random expander codes at the heart of our construction are not locally 
testable [1]. Interestingly, for our construction, better expansion translates into greater robustness. 
This should be seen in contrast to results in Ref. where good expansion (admittedly not 

of the same graph) appears to hurt the local testability of the quantum codes. We also remark 
that in the very recent result of |7] , quantum codes are constructed by applying Hi to classical LT 
codes, which leads to an alternative form of robustness where errors with small syndrome weight 
correspond to highly entangled states. 

The remainder of this extended abstract is organized as follows: Section describes the code 
construction with its expanding properties and states the main result. Section [^describes the basic 
decoding algorithm. Section gives the main points of the analysis of the algorithm, and states 
the robustness result. Sectiongives some concluding comments. The proofs of the more technical 
lemmas are relegated to appendices. 


2 Expansion, code construction and main result 

Let G = {AU B,S) he a biregular bipartite graph with left (resp. right) degree equal to (resp. 
A-b)- Let 1^41 = UA, \B\ = ub, and suppose ub < so that Aa < A^. We shall write a b 
(or more concisely, a ^ b when G is clear from context) to mean that the vertices a and b are 
adjacent in the graph G. If S' is a subset of vertices, denote by r(S) the set of all neighbors of 
vertices of S. Let us say that G is ( 7 ^, 5a)-Ig ft-expanding, for some constants 7 ^, 5a > 0, if for any 
subset S A with |5| < ^AnA we have |r(S)| > (1 — 5 a)Aa|S|. Similarly, we shall say that G is 
{iB-, 5b)- right-expanding if for any subset S ^ B with |S| < 'Jbbb we have |r(S)| > (1 — (5b)Ab|S|. 
Finally we shall say that G is {ja, 5a, 73 ,5b) left-right expanding (or simply expanding) is it is 
both ( 7 ^, (lA)-left-expanding and ( 7 s, (5B)-right-expanding. 

To any bipartite graph G we may associate the ns x ua matrix H, whose rows are indexed by the 
vertices of B, whose columns are indexed by the vertices of A, and such that Hij = 1 if i and j are 
adjacent in G and Hij = 0 otherwise. A binary linear code Cq can be thus dehned as the set of 
vectors x of such that Hx^ = 0, i.e. Cq is the code with parity-check matrix H. Conversely, 
any code C has a parity-check matrix H and the binary matrix H can in turn be viewed as an 
incidence relation between its rows and columns, i.e. a bipartite graph G: such a graph is usually 
called “the” Tanner graph or the factor graph of the code C. A code has several factor graphs, 
because it has several parity-check matrices, but it is usually clear which one we are talking about. 

Expansion and classical error-correction are connected through the following result of Sipser and 
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Spielman HZ]. 

Theorem 1. Let G = {AuB,£) he a (A^, /S.B)-hiregular ( 7 ^, 5a)- left-expanding graph. Letting 
and Ab be fixed and allowing ua to grow, there exists a decoding algorithm for the associated code 
Cg that runs in time linear in the eode length ha = \A\, and that, under the eondition 8a < 1/4, 
corrects any pattern of at most ^jAn-A errors. 

Our objective is to derive a quantum analogue of Theorem The codes we shall work with are the 
codes of m, whose construction we briefly recall. As mentioned in the introduction, a CSS code 
is defined by two classical codes Cx and Cz- in our case, both these classical codes are constructed 
from a fixed bipartite graph G = {AUB,£). Let us describe Cx and Cz through their factor graphs 
Qx and Qz- 

The bipartite graph Qx has left set of vertices U B^, and its right set of vertices is A x B. The 
bipartite graph Qz has the same left vertices but its set of right vertices is B x A. We will find 
it convenient to denote vertices of Qx and Qz by pairs of letters, omitting parentheses to lighten 
notation, and to use Greek letters for right vertices of Qx, Latin letters for right vertices of Qz, and 
denote elements of by a Greek letter followed by a Latin letter, and elements of B^ by a Latin 
letter followed by a Greek letter. Typical elements of AP‘,B^,A x B, and B x A will therefore be 
written respectively, aa, b/3, afi, and ba. The incidence structure of the two graphs is defined as 
follows: 

• In Qx ■ the set of neighbors of left vertex aa G is defined as 

r(Q;o) = {afi ^ A X B,a ~g fi}- 
The set of neighbors of left vertex 6/3 G B^ is defined as 

r(6/3) = {a/3 £ A X B, a ~g 6 }. 

• In the set of neighbors of left vertex aa G is defined as 

r(a;a) = {6a £ B x A, a ~g b}- 
The set of neighbors of left vertex 6/3 G B^ is defined as 

r(6/3) = {6a £ B X A,a ~g /3}- 

The corresponding parity-check matrices Hx and Hz of Cx and Cz may therefore be written in 
concise form as: 


Hx = , (1) 

Hz = {H ® Ing ® H^) . (2) 

It is not difficult to check [18] that the rows of Hx are orthogonal to the rows of Hz, so that 
(Cx,Cz) makes up a valid CSS quantum code. We also see that the parity-check matrices Hx and 
Hz are low density, with constant row weight A^ -|- A^. It is proved furthermore in [18] that the 
parameters of the quantum code Qg = (Cx,Cz) are 

lln = -h n^, k > (ua - nnf, min((i, d'^)]] (3) 
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where d denotes the minimum distance of the classical code Cg, i-e. the minimum number of 
columns of the ub x nA matrix H that sum to zero, and stands for the associated transpose 
minimum distance, which is the minimum number of rows of H that sum to zero. This is with 
the convention that each of these respective minima is set to equal oo if there are no subsets of 
columns, or rows respectively, that sum to zero. The transpose minimum distance associated to a 
parity-check matrix of a classical linear code is not a standard parameter because it is usually oo, 
but this is not necessarily the case for a matrix associated to a typical bipartite biregular graph, 
and it cannot be totally overlooked. 

Crucially, if one fixes the degrees Aa,^b with A a < and one takes an infinite family of 
graphs G = {AU B,£) with increasing number of vertices nA,nB, the rate of the quantum code 
Qg is bounded from below by a non-zero constant and its typical minimum distance scales as the 
square-root of its length n. 

We remark that the construction of [13] is somewhat more general, using two base graphs Gi and G 2 
rather than a single graph G. The above construction corresponds to the case when G = Gi = G 2 
and we choose to restrict ourselves to this setting for ease of description and notation. 

We can now state our main result: 

Theorem 2. Let G = {A\JB, £) be a (A^^, AB)-biregular ( 7 ^, 6 a , 7_b, 6b) -left-right-expanding graph. 
Assume the conditions 6 a < 1/6 and 6 b < 1/6. Letting A^ and Ab be fixed and allowing nA,nB 
to grow, there exists a decoding algorithm for the associated quantum code Qg that runs in time 
linear in the code length n = n\ + and that decodes any quantum error pattern of weight less 
than ^ 

Wo = ilAnA, Ibwb) • (4) 


Comments. Theorem implies in particular that the quantum code Qg must have a minimum 
distance proportional to min ( 7 ^ny 4 , 751 x 5 ), and that the decoding algorithm corrects a number of 
adversarial errors equal to a constant fraction of the minimum distance. We stress that we need 
both left and right expansion from the graph G, and that for these values of 6a and 6b, there are no 
known constructions that achieve it. However, the graph G with the required expanding properties 
may be obtained by random choice with classical probabilistic arguments, as developed in |3| or 

m- 


3 Description of the decoding algorithm 


We first need to describe the decoding problem precisely. An error pattern is a set of coordinates 
Ex C {1,... ,n} on which there is an X Pauli error, together with a set of coordinates Ez on 
which there is a Z Pauli error. There may be coordinates in Ex £iEz on which an X error and a Z 
error occur simultaneously, which is equivalent to a T Pauli error. The error pattern can therefore 
be given in equivalent form as the couple of binary vectors e = {ex,ez) where ex and ez have 
supports Ex and Ez respectively. The input to the decoder is the syndrome a{e) of the error e. It 
is made up of the classical syndromes of the vectors ex and ez for the matrices Hx and Hz, i.e. 
cj(e) = {crxiex),o'z{ez)), with ax{ex) = Hxe^ and crz(ez) = Hze^. We will say that the error 
pattern e = {ex, ez) is correctly decoded if on input a{e) the decoder outputs {ex + /x, ez + fz). 
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with fx and fz in the row space of Hz and the row space of Hx respectively. Our purpose is to 
exhibit a decoder that will always correctly decode every error e of weight smaller than the quantity 
given in Theorem]^ The weight w{e) of the error e = {ex,ez) is the number of coordinates i in 
which at least one of the two vectors ex Hz is non-zero. 

Before describing the decoder, let us fix some additional notation. The rows of the matrix Hz are 
called Z-generators (for the purpose of this paper’s internal notational coherence: note that this 
terminology is not standard in the stabilizer formalism). Recall that a row of Hz is indexed by 
an element 6 a of R x A. We choose to identify the corresponding generator with its support in 
A'^ U and denote it by gha, so that we have: 

gba ■= {aa : a ~ 6} U {b/3 : /3 ~ a} . (5) 

Similarly the X-generators are denoted by gag and we write 

9ag ■= {a/3 : a ~ /3} U {ab : 6 ~ a} . ( 6 ) 

The decoding algorithm treats X and Z errors independently. This means that the decoder applies 
two separate decoding procedures, one consisting of outputting ex + fx from crx{ex), and the 
other outputting ez + fz from aziez)- 

We remark that it is enough to prove Theorem]^ for error patterns of the form {ex, 0) and of the 
form (0, ez)- We will therefore describe the decoding algorithm for errors of type {ex,0), the other 
case being symmetric. Note that the weight of the error pattern w{ex,0) is simply the Hamming 
weight of the vector ex, and we denote it by \ex\- 

The decoding algorithm for ex works by going through the generators gba and for each one, by 
looking at whether flipping any pattern of bits strictly decreases the weight of the syndrome. More 
precisely: the algorithm takes as input a syndrome vector s G It then looks for a vector 

ei G such that: 

• the support of ei is included in some generator gba, 

• Si = s -|- o'x(ei), with |si| < |s|, 

• is maximal, subject to the above two conditions. 

If the algorithm cannot find a vector ei satisfying the first two conditions, it outputs a decoding 
failure. After i decoding iterations, the algorithm is in possession of a syndrome vector Sj G 
together with a vector ei -|- 62 + • • • + e* G such that 

• the support of every e^, j = 1 ... i, is included in some generator, 

• Sj = s -I- ax{ei -he,) and |si| < |s|. 

The {i I)-th decoding iteration consists in finding a vector ej+i G such that 

• the support of Cj+i is included in some generator gba. 
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• Sj+i = Sj + crx(e*+i), with |si+i| < |si|, 

• is maximal, subject to the above two conditions. 

The algorithm proceeds until it reaches some iteration i after which it cannot find any generator 
Qha which enables it to decrease the weight of s*. If |sj| / 0, then it outputs a decoding failure. 
Otherwise we have Sj = 0, and the algorithm outputs the vector eo := ei + 62 + • • • + e*. 

We shall prove that if s = crxi^x) with \ex\ sufficiently small, then the decoding algorithm never 
outputs a failure and its output eo satisfies eo + ex £ equivalently eo + ex is a sum of 
Z-generators, meaning the algorithm has correctly decoded the error pattern ex- 

By carefully updating the list of generators to be re-examined after every iteration, we obtain in a 
classical way an algorithm that runs in time linear in the number n of qubits, in the uniform cost 
model. 


4 Analysis of the decoding algorithm 

We first recall some tools and results from m- 

Consider the graph G = {Au B,S) and a subset of vertices S. Let us define the set r„(5) of 
unique neighbors of S, that is, the set of vertices v such that v has degree 1 in the graph S U r(S') 
induced by S. The complement of rtj(S') in r(5) will be called the set of multiple neighbors of S 
and denoted rm(*S'). Provided that the expansion in the graph G is large enough, then the graph 
displays unique-neighbor expansion, in the following sense: 

Lemma 3. Let G = {AU B,S) be a {'j a, h a)- left-expanding graph with 6 a < 1/2. Then, for 
any subset Sa ^ A with |5x| < '^a'^a, we have: 

ir^Sx)! > (l-2dA)AA|Sx|. (7) 

Proof. Consider the bipartite graph induced by Sa and r(5x). The number of edges incident to Sa 
should coincide with the number of edges incident to r(S'x), that is = |ru(S'y 4 )| + A|rm(*S'^)| 

where A > 2 is the average right degree on rm(*S'x). This implies that |rm(5'A)| < ^(AxIS'xl — 
|rM(5yi)|. Moreover, the expansion in G requires that |ru(5x)| > (1 — 5 A)Ax|5'y4| — |rm(*S'A)|. 
Combining both inequalities gives the unique-neighbor expansion. □ 

We also recall from m 

Proposition 4. If G is {jATd a)- left-expanding for 6 a < 1/2, then the minimum distance d of the 
associated classical code Cc is at least 'jAnA- 

From Proposition we immediately obtain: 

Corollary 5. If the graph G = {AllB,£) is a-, d a, -left-right-expanding with 6a, 6b < 1/2, 

then the minimum distance of the associated quantum code Qo is at least min( 7 xnyi, 
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Proof. Proposition implies that the minimum distance of the associated classical code is at least 
'jAnA, and, by inverting the roles of A and B, that the transpose minimum distance is at least 
'ysns- The result follows from (|^. □ 

The crucial part of the original decoding strategy of Sipser and Spielman consists in showing that, 
for errors of sufficiently small weight, there must exist (at least) one critical variable vertex of A 
in error that has many unique neighbors. Flipping the bit associated with this vertex decreases 
the syndrome value and this eventually leads to showing that one can always decode correctly by 
flipping bits that decrease the syndrome weight, provided the initial error is sufficiently small. In 
the present setting, we need a corresponding notion of criticality for generators. It will build upon 
the unique neighbor idea with a number of twists. 

Definition 6. Let E CZ A'^ U he a subset of vertices that can he thought of as an error pattern. 
Let us say that a generator g^a = {ota : a ~ 6 } U {6/3 : /3 ~ a} zs critical (with respect to E) if it 
can be partitioned as follows: 

Qba = Xa U Xa U Xa U Xft U Xfe U (8) 

where 

• Xq U Xq U Xa is a partition of g^a H and Xb U Xfe U Xb is a partition of gba H B^, 

• Xq U Xft C Fi, Xa n Fi = 0 and Xb H Fi = 0, 

• every vertex in r(X(i) n r(xb) has exactly two neighbors in E, 

• every vertex in r(xa) n r(xb) has no neighbor in E, 

• every vertex in r(xa) n r(xb) and every vertex in r(X(i) n r(xb) has exactly one neighbor in 

E, 

• XaUXb^^, \xa\ < 25 bAb,\ xb\ < 

Note that Definition [^implies in particular that the syndrome vector has value 0 in the coordinates 
indexed by r(xa) nr(xb) and r(xa) nr(xb), and has value 1 in the coordinates indexed by r(xa) H 
r(xb) and in r(xa) n r(xb). See Fig. [^for details. 

Lemma below asserts that a critical generator with respect to an error pattern always exists, 
whenever the error pattern E has sufficiently small weight (cardinality), and Lemma claims that 
its is always possible to modify the pattern of qubits inside a critical generator in a way as to 
simultaneously decrease the syndrome weight and not increase the error weight. Recall that we are 
decoding an X-error ex that we simply denote e from now on. 

Lemma 7. If the weight of the error e satisfies 0 < |e| < min( 7 ^nA, 75 ^, 5 ), then there exists a 

critical generator gba with respect to the support E of e. 

The proof of Lemma is given in Appendix 

If e is a vector of F 2 , let us call its reduced weight, denoted WR{e), the smallest Hamming weight 
of an element of the coset e + C^, i.e. the set of vectors of the form e + /, where / is a sum of 
Z-generators. 


Lemma 8. If the reduced weight of the error e satisfies 0 < wji{e) < m.in{'-fAnA,'yBnB), then there 
exists a critical generator together with a vector ei whose support is included in the generator, such 
that |c7x(e)| — |o'x(e + ei)| > |ei|/3. In words, flipping the k bits of the support of e\ decreases the 
syndrome weight by at least fc/3. Moreover, wji{e + ei) < wr^c). 

The proof of Lemma is given in Appendix 

The decoding algorithm is analyzed in two steps. First we show that it never outputs a decoding 
failure, i.e. always decreases the syndrome weight to zero, and secondly we prove the output error 
vector Co is equivalent to the original vector e modulo the space of Z-generators. 

Lemma implies the robustness of the quantum code. 

Corollary 9 (Robustness). Any error e with reduced weight WR{e) < min( 7 ^nA, 7 _Bns) has a 
syndrome with weight bounded from below as |cr(e)| > ^wji{e). 

Proof. Without loss of generality, suppose e is a representative of e + with minimum weight. 
As long as the weight of the syndrome is positive. Lemma [^guarantees the existence of a generator 
and a vector ei whose support is included in the generator such that |crj!s:(e)| — |crx(e + ei)| > |ei|/3. 
Moreover, since wji{e+ei) < |e|, we can apply Lemmaj^again to e+ei and iterate, say i times, until 
the syndrome weight reaches 0. After i iterations, we obtain that the syndrome of e+ei+e 2 + - • -+6* 
is 0, hence that 

kx(e)| > ^(|ei| + |e 2 | H-|ei|) > ^|ei H- ei\, (9) 

and that WR{e + ei + • • • + Cj) < |e| < min( 7 ^nA, 7 _BnB). This last fact implies by Corollary [^ 
that e + ei + • • • + Cj is in C^, i.e. that e is equal to ei + • • • + modulo a sum of Z-generators. 
Inequality (|^ proves therefore that \ax{e)\ > ^wji{e). □ 

Unfortunately, the decoding algorithm is not guaranteed to follow the good decoding path exhibited 
by Lemma Indeed, the decoding algorithm simply tries to optimize the weight of the syndrome, 
but the error weight might increase in the process. Nevertheless, we have: 

Lemma 10. If the Hamming weight of the initial error e is less than wq, then the decoding algorithm 
never outputs a decoding failure and always correctly decodes e. 


Proof. The decoding algorithm chooses a sequence of vectors ei, 62 ,..., such that 

Wxie)\, \cFx{e + ei)|,..., jaxie + ei H-h ei)\,... 

is a decreasing sequence. Set eo = e, ei = e + ei,..., e* = ej_i + e*. Now whenever 

\£i\ < min(7AnA,7BnB), 


Lemma (HI ensures that 


By hypothesis (10) holds for i 


Wx{£i)\ - \(7x{£i+l)\ ^ 1 
|Ci+l| 3 

0. Suppose it holds for 0,1,..., i, we then have 

|ei| < 3(|(Tx(eo)| - kx(ei)l) 

... < ... 


|ei+i| < 3(|CTx(£i)| - |o-x(ei+i)|) 


( 10 ) 

( 11 ) 
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By summing all these inequalities we obtain |eiI+ ••- + 164+11 < 3(|crx(eo)| ~ |c’'x(£j+i)|) < 3|(Jx(e)|. 
Writing |6j+i| < |e| + |ei| + -•- + 164+11 we get |e 4 +i| < |e|+3|(Tx(e)| < |e|+3As|e| since the syndrome 
iTx(e) has weight at most Ab|6 |. From the hypothesis on the Hamming weight of 6 we get that (10) 
holds for z + 1, and inductively until the eventual output of the algorithm 6o = ei + - - - + 6j, such 
that £j = e+6o has zero syndrome. Condition (10) translates to |e + 6o| < mmljAnA, 'ysns) which 
ensures by Corollary that 6 + eo ^ C^, which means exactly that eo is a valid representative of 
the initial error e. □ 


5 Concluding comments and questions 


We have exhibited a linear-time decoding algorithm that corrects up to adversarial quantum 

errors over n qubits. While this is the largest such asymptotic quantity to date, one would hope 
to break this barrier and eventually achieve correction of H(n) errors. If one were to do this with 
quantum LDPC codes, this would imply obtaining the elusive proof of existence of low-density 
codes with a minimum distance scaling linearly in the number of qubits. 

Kovalev and Pryadko have shown |15j that the codes of [18] have the potential to correct number 
of random depolarizing errors that scales linearly in n, with a vanishing probability of decoding 
error. This is without decoding complexity limitations however, and a natural question is whether 
the ideas of the present paper can extend to decoding (n) random errors in linear or quasi-linear 
time. 

We have worked to achieve the smallest possible value of Sat^b in Theorem]^ i.e. the smallest 
possible expansion coefficient for the base graph G. Can the bound SA,dB < 1/6 be decreased 
further ? A somewhat related question is whether the left-right expanding base graph G can be 
obtained constructively, rather than by random choice ? 
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Appendices 



<■ 


X, 


> <■ 


X, 


> < - 


Figure 1: A critical generator gi,a- The generator g^a is partitioned into 6 sets: Xa and Xb contain 
the errors we want to correct, Xa and yb do not contain any error, and Xa and Xb have a small 
relative size. When flipping the qubits corresponding to Xa and Xb (or the qubits in Xa and Xb), 
the syndromes in ^ab see their weight decrease from 1 to 0 while the syndromes in the sets 

Sab and S^b remain unchanged; hnally either the syndrome weights in 5a U 5^, or in Sa U Si can 
increase by 1 while the others stay unchanged. 


A Proof of Lemma 0 


First, we need the following extra notation. Let 5 be a subset of A in the graph G = {AuB,£): For 
a vertex a G 5, we denote as follows, respectively, the set of unique neighbors of o in the subgraph 
of G induced by 5 U r(5), and the set of multiple neighbors of a in the same induced graph: 

rf(a) :=r(a)nF„(5), r^(a) := F(a) n r™(5). (12) 

Denote now by 7 ^ 0 the support of the X error e to be corrected and dehne Ea ’■= E A? 
(resp. Eb := E r\ B^) the set of error vertices in A^ (resp. in B^). 

Suppose hrst Eb = 0. Consider the projection Ej^ of Ea on the second coordinate: 

E\ := {a G a : Ax {a} H Ea 7 ^ 0} . (13) 


Clearly \EW < \Ea\ < w{e) < XAnA- Lemma then implies that |ra(F^^)| > (1 — 25 a)Aa|F^^|. 
Therefore, using the notation (12), there exists a G E\ such that |ra'^(a)| > (1 — 25a)Xa- Pick 
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some aa G Ea and b ^ a. Then the generator gi,a can be partitioned into four components 
9ba = Xa U Xa U Xb U Xb which satisfy the requirements of the lemma. Specifically, we define 
Xa = 9ba n Ea, Xa as the complement of Xa in gba H A^, 

Xb = {b/3egbanB^ : /3 G r5(a)} 

and Xb as the complement of Xb in la words, Xa and Xa partition the part of the generator 

gba into error and error-free components, and Xb and Xb partition the part into, respectively, the 
set whose second coordinate is a unique neighbor and the set whose second coordinate is a multiple 
neighbor with respect to E\. 


The proof is identical with the roles of A and B interchanged if Ea = 0- 

Let us turn our attention to the remaining case where both Ea 7 ^ 0 and Eb 0- We define E\ as 
above and choose a coordinate a E E\ such that Tu^{a) > (1 — 25^)A^. From now on, we use 
the shorthand r„(a) for this set, leaving E'j^ implicit. Next, we define: 

EB,a ■= {bp G Eb^ P G r„(a)}. 

If EB,a = 0, we proceed as in the case where Eb = 0. Otherwise, by assumption on the weight of 
the error, namely |e| < xsnB, Lemmaapplies again, this time to the set Ej^ ^ corresponding to 
the projection of EB,a on its first coorcimate, 


rn(^) := we define: 


* 

{b : 

3/3 G I 

P{a),bp 

G Eb} ■ 

(14) 

Eh, a 

such 

that 


> (1 — 2(3s)As. 

Using the shorthand 

Xa ■ = 

{aa 

G Ea 

cx E F^/ 

m 

(15) 

Xa ■ = 

{aa 

^Ea 

oi E 

m 

(16) 

Xb ■■ = 

{bp 

G Eb 

/3er, 

(«)} 

(17) 

Xb ■■ = 

{bp 

^Eb 

^ef. 

(«)}■ 

(18) 


The construction ensures that every element a/3 that belongs both to the neighborhood of Xa U Xa 
and of Xb U Xb is such that a is a unique neighbor of b and p is simultaneously a unique neighbor 
of a, unique neighborhood being understood with respect to E. Moreover Xb / 0 by construction 
of the set E{^ a / 0- 


B Proof of Lemma [8] 


We consider the generator gba promised by Lemma Without loss of generality, we may suppose 
that the error vector e is in reduced form, i.e. it is the vector of lowest Hamming weight in the 
coset e -|- . The consequence is that the number of qubit vertices in error | E n gba \ within the 

generator gba is not more than (A^ -|- Ab)/ 2, otherwise replacing e by e -|- e', with e' the vector 
having gba for support, would yield a vector with strictly smaller Hamming weight within the coset 
e + C^. In particular we have: 

\Xa U Xb\ < + As). (19) 
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Let us introduce the following reduced variables: 


X := 


y ■■= 


\Xa\ 

Ab’ 

kbi 


\Xa\ 

X := 1 — X — z, 

(20) 

><\<\ 

II 

y --l-y -t. 

(21) 


By assumption, x + y > 0. 

The expansion condition Sa,Sb < g implies that z,t < 1/3 by the last condition of Definition 
However, zAb and xA^ are integers, this implies 3zAb < 3Ab — 1 and 3tA^ < 3Ayi — 1, which 
in turn give: 


o<.a 


1 

3Ab ’ 


0<t<l 


1 


( 22 ) 


We are looking for a vector ei, whose support is included in gf,a, such that 


i. the syndrome of e + ei has strictly smaller weight than the syndrome of e, and the difference 
|<7x(e)| — |o'x(e + ei)| is at least |ei|/3, 

ii. the reduced weight WR{e + ei) is at most equal to the Hamming weight |e| of e. 

We will consider four cases: 


1. ifx + y<2/3, then the vector ei is chosen to have support Xa U xj,, 

2. if X < X and y < y, then the vector ei is chosen to have support Xq U x^, 

3. if the above two hypotheses do not hold and if x > x and y > y, then the vector ei is chosen 
to have support either Xa U Xb, or its complement in g^a, 

4. in the remaining cases, we show that either Xq U Xb or Xa U Xb, or x^ U xj, U Xa U Xb is an 
adequate choice for the support of ei. 

Let us introduce the partition of T{gba) induced by the partition ([^: 

r(yba) = SaUSbUSaUSlU Sab U S^l U Sab U U 5 

with 

Sa =T{Xa)nT{xb) Sb =T{xb)nr{Xa) Sa =T{Xa)nT{xb) S-, =T{xb)nT{Xa) 

Sab = r(xa) n r(xfe) s^i =r(xa)nr(xb) Sab =r(xa)nr(xb) s^j, =r(xa)nr(x6) 

5 =r(xa)nr(xb) 

as represented on Figure 

Let us denote by d the decrease of the syndrome weight when we flip Xq U Xb (i.e. choose ei to 
have support Xa U Xb). Similarly we denote by d the decrease of the syndrome weight when we flip 
either Xa U Xb or Xa U x;, U Xa U Xb- 
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These quantities satisfy: 


(?/(AaAb) > xy + xy - xt - yz (23) 

d/{AA^B) > xy + xy — xt — yz. (24) 

This is because when we flip the bit values on the set U Xfe, then the value of the syndrome is 
changed from 1 to 0 on the support and S^b, remains at 0 on the supports Sab and may 
possibly change from 0 to 1 on the supports SaD Sb, and remains unchanged in the other regions. 

Similarly, when we flip the bit values on the set Xa U Xb, the syndrome is flipped from 1 to 0 on the 

support Sab and Sab and possibly from 0 to 1 on 5^ U S^ while the other regions remain unchanged. 

We now address the four cases: 


1. Suppose X + y < 2 / 3 . Then (23) can be rewritten as: 

d 


AaAb 


> xy + xy — xt — yz 

> x{l — y — t) + y(l — X — z) — xt — yz 

> X + y — 2xy — 2{tx + yz) 




by applying Eq.(22) to z and t. The first term is nonnegative for x + y < 2/3, so that we 
obtain: 

^ 2 , , 1 , 

O > -\Xa U Xb\ > -\Xa U Xb\. 

Furthermore, when ei has support Xa U Xb, the support of ei is included in the support of e, 
hence |e -(- ei| < |e| and WR{e + ei) < |e|. 


2. If a: < X and y < y, then y = l — t — y<y gives y > 2 — 21 hence, applying (22), 


y-t> 


Similarly we have x — z > and (23) gives 


2Aa' 


d > AaAb [x{y -t) + y{x - z)] 

> - [xAb + yAA] 

1 , 

> - |Xa U Xb\ . 

The set Xa U Xb is again an adequate choice for the support of ei, and tc/j(e + ei) < |e| as 
before. 

3. If X > X and y > y, then by an argument symmetrical to the preceding case, exchanging x 
with X and y with y, we get: 

9 > ^ l^a U Xb\ ■ 
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We remark that we must have \xa U xi\ > 0, otherwise z,t < 1/3 imply that \xa U Xh\ > 
(A /1 +As)/2, which contradicts our assumption (19). Therefore, choosing ei to have support 
XaUxb gives a strict decrease in the syndrome weight that is at least |ei|/3. Now in this case, 
the Hamming weight |e + ei| is larger than |e|. However we need to consider its reduced weight. 
A vector equivalent modulo to e + ei is e + e/ with the support of e'l being XaUxbUXa^Xb- 
The Hamming weight of e + e'l is 

|e + e/I < |e| - (x + y) + {z + t). 

But since x + y > 2/3 and z, t, < 1/3, we have |e + e/| < |e|, so that wji{e + ei) < |e|. 

4. Finally, in the remaining case we consider 9 + 9 to show that one of the two quantities, d or 
8 is sufficiently large. Adding (23) and (24) yields 

{8 + 8) > 2xy + 2xy — t(l — z) — z(l — t). 


We now observe that the negation of the condition 

{x < X and y < y) or {x > x and y > y) 
implies 2xy + 2xy > (x + x){y + y) = {1 — z){l — t), from which we get: 

1 


AaAb 


(5 + 9) > (1 — z){l — t) — t(l — z) — z(l — t). 


The formal equality 

(1 — z){l — t) — z{\ — z) — t{l — t) = 3 


— z 


-t] - 


1 

3 


yields therefore, when combined with (22), 

1 , „ 1 


Hence, after rearranging, 


1 


(9 + 5)>- l + — 1 + — -- 


Aa 


1 


As 


1 A^ + As 


max {a, 9} > -((9 + 5) > - ^ 

if max|5, 5} = 5, then we choose ei to have support XaUx;,, and we get 8 > ^|ei| from 
condition ( [I^ . We have WR{e + ei) < |e| as in cases 1. and 2. 

if max {5, 5} = 5, then 

— either |xa U Xfe| < ^ which case we set ei to have support Xa U Xb and have 

8 > jleil and wr^c + ei) < |e| by the same argument as in case 3, 

— or jxa UXftI > ^ which case we set ei to have support XaUXbUXa^Xb so as 

to again have 5 > g|ei|. The Hamming weight |e + ei| is at most \e\ — {x+y) + {z+t), 
which is less than |e| as in case 3, so that again WR{e + ei) < |e|. 
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